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We have designed and characterized micro-electro-mechanical systems (MEMS) for applications 
at low temperatures. The mechanical resonators were fabricated using a surface micromachining 
process. The devices consist of a pair of parallel plates with a well-defined gap. The top plate 
can be actuated for shear motion relative to the bottom fixed plate. Details on the operation and 
fabrication of the devices are discussed. A full characterization in air and at room temperature was 
conducted as a function of pressure, from 10 mTorr to 760 Torr, ranging from a highly rarefied gas to 
a hydrodynamic regime. Finally, we demonstrate the operation of our resonator at low temperatures 
immersed in superfluid 4 He and in the normal fluid of 3 He down to 50 mK. These MEMS oscillators 
show potential for use in a wide range of low temperature experiments, in particular, to probe novel 
phenomena in quantum fluids. 

PACS numbers: 07.20.Mc, 85.85.+j, 67.10.Jn 



Mechanical oscillators such as torsional oscillators, vi- 
brating wires, and quartz tuning forks have been used 
extensively at low temperatures for the study of quan- 
tum fluids. They allow the direct probe of properties of 
liquid helium and have provided irrefutable evidence of 
superfluidity These devices can be built in- house or 
are commercially available, providing a fair amount of 
flexibility in designing experiments. Modern silicon tech- 
nology allows on-demand fabrication of micro-electro- 
mechanical systems (MEMS) with a plethora of function- 
alities as sensors and actuators. The capability to create 
devices and tailor its mechanical and electrical proper- 
ties provides a rich arena in which novel experiments can 
be designed to explore new physics in quantum fluids at 
reduced length scales. 

In recent years, micro/nano-mechanical structures 
have been incorporated in the low temperature study to 
reach the quantum mechanical ground state of mechani- 
cal resonators 0, Q or to create devices similar to vibrat- 
ing wires with the intention of using them for the study 
of liquid helium 04^] • Predictions of a novel superfluid 
state in confined 3 He Q have also stimulated the devel- 
opment of experimental probes capable of confining and 
probing the liquid's properties at the nanometer scale [1- 
1 101 ] . Our motivation to develop the devices described in 
this paper stems from our interest in investigating novel 
phenomena occurring in superfluid 3 He films. We envi- 
sioned a device that can form well-defined films of liquid 
and can simultaneously probe their properties with a high 
resolution. Our devices meet these design criteria. 

In this paper we describe the design and characteriza- 
tion of the MEMS oscillators. Details on the design, fab- 
rication, and preparation are described in section I. The 
morphology of the surface is an important quality which 
could have an influence on the interaction of the device 
with its surrounding medium. We characterized the sur- 
face through atomic force microscopy (AFM) measure- 
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FIG. 1: 

CAD image of a typical HI device. No substrate is shown in 
this figure. The springs also serve as electrical connections 
from the center plate to the external boding pads on the 
chip. 



ments. The results from these investigations are shown 
in section II. The general properties of the devices and 
the detection scheme are described in section III. The 
device was studied at room and low temperatures. At 
room temperatures, the resonant properties were stud- 
ied through a wide range of pressure, from 10 mTorr to 
1 atm. Finally, we demonstrate the use of a device below 
1 K while submerged in the superfluid state of liquid He 
and in the normal state of liquid 3 He down to 50 mK. 
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TABLE I: Thickness of the different layers used in Poly- 



MUMPs 

Layer name Material Thickness (jim) 

Nitride Silicon nitride 0.6 

PolyO Polysilicon 0.5 

First Oxide PSG 2 

Polyl Polysilicon 2 

Second Oxide PSG 0.75 

Poly 2 Polysilicon 1.5 

Metal Gold 0.75 



I. MEMS DEVICE FABRICATION AND SETUP 

The geometry of the device is that of the so-called 
comb-drive actuator. It is similar to the one used in mod- 
ern MEMS accelerometers G3 • The device consists 
of a movable center plate («200x200 /im 2 ) suspended by 
four serpentine springs. The center plate is maintained 
at a fixed distance above the substrate, thus creating a 
uniform gap. The two devices used throughout this work 
have the gaps of 1.25 and 0.75 /im, and are referred to 
as HI and H2, respectively. A CAD image of a typical 
HI device is shown in Fig. [T] The center plate can be set 
in lateral motion by electrostatic interaction between the 
inter-digitated comb-like electrodes attached to its sides 
and the fixed electrodes anchored to the substrate. For 
example, if a dc voltage is applied across one pair of elec- 
trodes, the center plate will be pulled towards the fixed 
electrode along the x-direction (see Fig.QJ. 

We employed a shared wafer surface micromachining 
process called PolyMUMPs. The chips were fabricated 
by a commercial foundry called MEMSCAP. A layout of 
the different lithographic levels to be patterned on the 
chip was done using a specialized CAD software called 
L-Edit. The designs were then sent to the company for 
fabrication. In this process, the following layers are de- 
posited on the substrate: three polysilicon layers (PolyO, 
Polyl, Poly2) used as the structural material, one sil- 
icon nitride layer (Nitride) used for electrical isolation 
from the substrate, two phosphosilicate glass (PSG) lay- 
ers (First and Second Oxide) as sacrificial layers, and a 
metal layer (METAL) used for improved electrical con- 
nectivity. The label, material, and thickness of each layer 
patterned on the wafer is listed in Table |TJ The pro- 
cess is carried out on 100 mm wafers heavily doped with 
phosphorus. The layers are deposited using low pressure 
chemical vapor deposition (LPCVD) and patterned by a 
combination of photolithography and reactive ion etching 
(PJE). PSG is used in the process to dope the different 
Poly layers with phosphorus by annealing the wafer at 
1050°C in argon during different steps of the process. 
The typical resistivity of a Poly 1 layer is 2.61 x 10~ 3 £1- 
cm. After depositing and patterning all the necessary 
layers to create the structures, the resulting chips are 
sent to the customer. More details can be found in the 




FIG. 2: 

SEM pictures of a typical device. An overall view of an H2 
is shown in the top left figure. The top right figure shows 
half of the same device. The bottom left figure shows a more 
detailed view of the spring of an HI device at its anchor 

point with the substrate. A close up look at the 
comb-electrode fingers is displayed in the bottom right 
figure. 

design handbook at the company's website [l3j . 

The die received from the foundry was diced into 
smaller square chips of ~ 2.5 mm side length using a 
diamond saw. During this procedure, the devices were 
still coated by a photoresist layer, and the sacrificial layer 
was still present. This prevents damage to the device be- 
cause the polysilicon structures are not yet movable. Af- 
ter dicing, the devices were "released" using a wet etching 
procedure. This removes the sacrificial PSG layers and 
frees the mechanically active parts. A standard process 
suggested by the foundry is followed for releasing [13j. 
After the release process, if the chips are allowed to dry 
by exposure to the environment, the small plates com- 
prising the devices might be pulled towards each other 
by capillary action and stick together. This phenomenon 
is called stiction. To avoid this problem the chips were 
dried using a CO2 critical point dryer. After drying, the 
chips were packaged in a custom designed 20-pin socket. 
Electrical connections were made from the bonding pads 
on the chip to the pads on the package by using a West- 
Bond ball wedge wire bonder with 0.025 mm thick gold 
wire. Scanning electron microscope (SEM) images of a 
typical released device are shown in Fig. [2j 



II. AFM SURFACE CHARACTERIZATION 

Experiments in thin films of 3 He have found a possible 
non-trivial relation between the roughness of the surface 
bounding the liquid and the scattering relaxation time 
of quasiparticles close to the surface [l4[. A theoretical 
description developed by Meyerovich et al. transforms 
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FIG. 3: 

Top: Topography of the surface of the Poly layer in a 
typical device. Bottom: Histogram of heights from the 
topography of the surface of a device. The solid line 
represents a fit to a Gaussian. 



FIG. 4: 

Top: Normalized autocorrelation function calculated from 
the topography of the Poly surface in a typical device. 
Bottom: ID normalized autocorrelation function averaged 
from different spots on the device. 



the problem of electrons propagating through a film with 
rough surfaces and a homogeneous bulk into that of a 
film with flat surfaces and a disordered bulk p^| - [l7| . This 
theoretical scheme was adapted to explain recent exper- 
iments in which an anomalous temperature dependence 
was observed in the transport properties of a thin film 
of liquid helium |18j |. This theory uses, as an input, in- 
formation about the surface's autocorrelation function, 
which is related to the momentum transition probability 
of quasiparticles scattering off the surface. It is, there- 
fore, imperative to have a detailed understanding of the 
topography of the polycrystalline plates in our devices. 
For this purpose, we have carried out AFM measure- 
ments on a sample chip containing the MEMS resonators. 

The chip was initially sprayed with nitrogen gas to 
blow off the top plate of all the devices present on the 
chip. A Nanoscope III A Scanning Probe Microscope was 
used to scan the exposed Poly surfaces that consti- 



tute the bottom plate of a typical device. Twelve differ- 
ent spots were probed with scans ofl0xl0,2x2, and 
0.4 x 0.4 /L*m 2 on different parts of the same device and on 
different devices contained in the same chip. The AFM 
tip was set in tapping mode with a resonance frequency 
of 69.16 kHz. Data for a scan size of 2 x 2 /im 2 are 
shown in Fig. [3] The data was flattened and plane fit- 
ted. A histogram of heights shows a symmetric Gaussian 
distribution about the mean plane, which is indicated as 
zero height. The width of the Gaussian is 20.8 nm, thus, 
the average feature height with respect to the mean plane 
is 10.4 nm. 

A normalized autocorrelation function was calculated 
from the surface topographies to investigate the corre- 
lation between grains. As shown in Fig. [4J the func- 
tion is sharply peaked at zero displacement. From the 
autocorrelation corresponding to the 2x2 fim 2 topog- 
raphy, a one-dimensional autocorrelation was calculated 
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TABLE II: Calculated device properties from layout geometry 



Device 


Gap, d 


Area, A s 


Mass, m 


Spring, k a 


Capacitance, Co 




w 


(/jmX(im) 


(xHT 10 kg) 


(x4 N/m) 


(xl(T 3 pF) 


Hlt.m 


1.25 


196 x 196 


3.45 


2.47 


8.41 


H2 


0.75 


178 x 178 


2.77 


2.47 


8.41 




FIG. 5: 

Circuit diagram for the capacitance bridge technique employed. An external ratio transformer is used as a tunable inductive 
voltage divider along with the two capacitors formed by the two sets of electrodes at the sides of the device. 



by averaging all the values on the perimeter of a circle of 
radius r centered at the origin of the contour plot (top 
panel of Fig. [4j . The resulting autocorrelation function is 
shown at the bottom of Fig. HI The autocorrelation goes 
from a positive to a negative value at point Pi, where 
r = 0.137 /im. This distance provides an estimate of the 
average grain lateral size. At the second point, P2, the 
autocorrelation goes back to a positive value. This point 
indicates the existence of another feature on the surface 
of average size 0.591 /im. Clusters of grains of approx- 
imately this size can be seen in the topography of the 
surface. In summary, from the height histogram, grains 
can be estimated to have a height of 10.4 nm. From the 
autocorrelation, the estimated width is 0.137 /im. Fur- 
thermore, an average cluster size can be estimated to be 
0.591 /im from the autocorrelation function. The last two 
distances establish two length scales: a short length scale 
at 137 nm and a long length scale at 591 nm. These can 
be verified by visual inspection of the topography (see 
Fig.©. 



III. DEVICE PROPERTIES AND 
CHARACTERIZATION 

Some important device specifications for both types of 
devices, HI and H2, are shown in TableUT] These include: 
the gap distance d, the area of the main plate ^4^, the 
total mass of the movable plate to, the spring constant 



k Sl and the capacitance cq. The total capacitance, Co, 
between a fixed electrode and its movable counterpart 
can be calculated from the expression [19( : 



Co 



Aee 



(1) 



where A is the total overlap area, df is the distance be- 
tween fingers, n is the total number of fingers, eo is the 
permitivity of free space, and e is the relative permitivity. 

The spring constant is determined by Young's modu- 
lus of polysilicon and the geometry of the spring. For a 
serpentine spring with N guided beams (N = 2 in our 
devices), the spring constant is given by [2fJ 



*(t) S ' 



AT' 



(2) 



where E is Young's modulus, w is the width, t is the 
height, and / is the length of the spring beam. The litho- 
graphic masks used in the process by MEMSCAP have a 
pixelation of 0.25 /im (l3| . If an error as high as 0.5 /im 
per lateral dimension is assumed, the total error in the 
spring constant can be as high as 38%. The mass of the 
device was calculated directly from the geometry of the 
devices. In this case the error in the lateral dimensions is 
negligible, and the only relevant source of the error is the 
thickness of the Poly layers (within ~ 10 nm) This 
introduces an error of less than 1% in the estimation of 
the mass. The resonance frequencies obtained from the 
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TABLE III: Resonance frequencies and measured Q-factors in high vacuum (~ 3 mTorr) 



Device 


j c (calculated) 
(Hz) 


js (simulated) 
(Hz) 


jm (measured) 
(Hz) 


Q — j actor (measured) 


HI 


29133.3 


24755.4 


23063.1 


6424.3 


Hl t 


29133.3 


24755.4 


22006.2 


23188.8 


Hl m 


29133.3 


24755.4 


21516.4 


105472.6 


H2 


32513.2 


28438.3 


26691.6 


5858.6 




FIG. 6: 

Vibrational modes for device HI. 



values of the calculated mass and the spring constant 
are listed in Table IIIII We used w = 3 ^m, t = 2 /im, 
I = 120 ^m, and E = 158 GPa. 



A. Detection Scheme 

The top plate is electrostatically coupled to the fixed 
electrodes through comb-like structured terminals. This 
allows the device to be actuated in shear motion, paral- 
lel to the substrate. The two fixed electrodes on oppo- 
site side form a pair of series capacitors with the center 
plate as a common electrode, which allows us to imple- 
ment a differential capacitance measurement by form- 
ing a bridge circuit with an external ratio transformer 
(TEGAM 1011 A). A schematic diagram of the measure- 
ment scheme is provided in Fig.[5j A carrier signal Vhf, at 
150 kHz, is provided by the internal oscillator of a Signal 
Recovery 7124 lock-in amplifier. The ratio transformer 
splits this carrier signal into opposite phases that are 
connected respectively to the two fixed electrodes. Vibra- 
tions of the top plate are induced by another ac voltage 
at a much lower frequency. The low frequency excitation 
(Agilent 33220A), Vf, is added to Vhf through a home- 
made summing amplifier. The frequency of Vif is swept 
through the resonance. As the device goes through reso- 



nance the off-balance signal is picked up by an Amptck 
A250 charge sensitive amplifier. The amplification fac- 
tor, a, of the A250 is given by the inverse of the feedback 
loop's capacitance c/ = 1 pF. The signal is later demod- 
ulated by two lock-in amplifiers connected in series. The 
first lock-in amplifier (Signal Recovery 7124) demodu- 
lates the carrier signal and the second lock-in (Signal 
Recovery 5210) detects the amplitude of the oscillations 
of the plate at the driving frequency (Vif) or its multi- 
ple. A dc bias voltage Vb was applied to one of the side 
electrodes of a device through a 2 kSl resistor. The total 
electrostatic force on the actuation side of the device can 
then be calculated and is given by 

F e (t) = ^(V b + V Q smutf (3) 
Idf 

= ^[(n 2 + -^)+2^osin^-^cos2^ . 

This force can be divided into three components: Fq, 
Fi, and F 2 . The first term, F oc V^ + V^/2, is the 
constant term in parenthesis on the second line of Eq. [3] 
The second term, F\ oc 2Vf,Vb sinwi, is in phase with the 
excitation and exists only when a DC bias is applied. 
v 2 

The last term, F2 oc cos 2wt, is at twice the applied 
frequency. F± can be written as 

Fy(t) = ^HVb smut = PVac{t), (4) 

where /3 = nteeoVb/df is often called the electro- 
mechanical transduction factor. 



B. Frequency Sweeps 

After a device was wirebonded and the circuit was set 
up, the resonance peaks were found through frequency 
sweeps while keeping the device in a closed chamber un- 
der vacuum, at a pressure P ~ 3 mTorr. The amplitude 
and resonance frequency were determined from fitting to 
a Lorentzian curve. To connect the mechanical and the 
electrical systems, the transduction factor in Eq.@]is used 
to couple Fi(t) and V ac (t) as well as i(t) and the velocity 
of the movable plate, v(t): 
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FIG. 7: 

Top: Amplitude and resonance frequency as a function of bias voltage for the shear mode of an HI device. Bottom: 
Amplitude and resonance frequency as a function of bias voltage for the pivot mode of an HI device. The lines represent a fit 
to a parabola. The resulting equations for the fit are displayed in each graph. 



The previous expression also establishes a direct rela- 
tion between the charge, q(t), and the displacement, x(t), 
through the transduction factor. A typical value for this 
factor is j3 » 1.8 x 10~ 9 Vb C/m. The in-phase (q x ) and 
out-of-phase (q v ) components of the charge detected can 
be found to be: 



q y {uj) 



(w 2 - w 2 ) 2 + cjAcj 2 

gojuAuj) 2 
(w 2 - w 2 ) 2 + (wAo;) 2 



(6) 



where go is the amplitude, ujq is the natural frequency of 
the resonator, and Aw is the width of the peak. In all 
of our devices tested, two resonance peaks were found. 
These correspond to different modes of oscillations. One 
of the peaks was found at ~ 14 kHz and the other at 
~ 23 kHz. The position of the peaks changes from de- 
vice to device due to the variations in the fabrication 
process as explained above. In order to identify the type 



of mode corresponding for each resonance frequency, nu- 
merical simulations based on finite element analysis were 
performed using COMSOL Multiphysics. The results for 
the different modes and their corresponding frequencies 
for HI are summarized in Fig. [5] For these simulations, 
a fixed boundary condition (no motion) was used at the 
anchor point of the spring (see Fig. [T]). Using the infor- 
mation from the simulations, the two modes found from 
the frequency sweeps were identified as the pivot mode 
about y-axis for the low frequency peak and the shear 
mode for the high frequency peak. The results of fre- 
quency sweeps, simulations, and analytical estimations of 
the resonance frequencies are summarized in Table HTT1 for 
the shear mode. Three different HI devices were stud- 
ied, indexed as HI, Hl t , and Hl m . The dimensions of 
these devices are the same, but the designs of the electri- 
cal connection from the bonding pads to the electrodes 
are different. For devices labelled with subscript "t", the 
connections were made "thicker" , while for the ones la- 
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belled with subscript "m", the electrodes were covered 
with a gold layer to improve conductivity. 



DC Bias Test 



In the absence of a dc bias, the second term in Eq. [3] 
(iq) would disappear and no peak would be detected 
when attempting to detect a signal at the same frequency 
as the driving force. However, it was found through these 
studies that even at Vb = there was still a resonance 
peak. This means that there is an intrinsic bias voltage 
between the center electrode and the driving electrode. 
As of now, we have not identified the source of this in- 
trinsic voltage. The detected resonance peak was found 
by scanning in "1/" mode in the second lock-in amplifier, 
which means the detected signal and the driving signal 
are at the same frequency. 

The bias voltage was changed from -10 to 10 V and the 
frequency was swept through resonance. The results for 
an HI device with metallic electrodes (Hl m ) are shown 
for the shear mode (top panel of Fig. [7]) and for the pivot 
mode (bottom panel of Fig. [7]). The amplitude increases 
linearly, as expected from Eq.[3j until it starts to saturate 
and to turn back down. This is a direct consequence 
of non-linearity starting to affect the peak due to high 
driving amplitudes. The resonance frequency follows a 
parabolic dependence on DC bias. The frequency of the 
shear mode decreases with the bias, while the frequency 
of the pivot mode is seen to increase with the bias. 

Since the frequency of the F2 component of the elec- 
trostatic force in Eq. [3] is twice of the excitation signal, 
the lock-in amplifier was set to detect vibrations at twice 
the input frequency ("2/" mode) and the frequency was 
swept around half the value of the original position of the 
peak. Using the amplitude of the peak in "2/" mode, it 
is possible to obtain a measure of the intrinsic bias by 
taking a ratio with the amplitude in "1/" mode with no 
bias voltage applied. The ratio of the amplitudes for an 
H2 device is found to be 



Af _ 2V b V _ 0.792 (mV) 
~~ V 2 /2 ~ 0.117 (mV) 



6.76, 



(7) 



from which we get \% = 0.425 V using Vb = 0.25 V. Sim- 
ilarly, using the shear mode peaks at "1/" and "2/", we 
found Vb = 0.482 V. These values are in good agreement 
with the results obtained from the parabolic fit shown in 
Fig. [7] Similar values are also obtained for all devices on 
different chips. 
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FIG. 8: 

Top: Q-factor and shear mode frequency of device Hl m as a 
function of pressure. Bottom: Q-factor and shear mode 
frequency of device H2 as a function of pressure 



IV. ROOM TEMPERATURE AND CRYOGENIC 
CHARACTERIZATION 

A. Resonance as a Function of Pressure at Room 
Temperature 

The influence of air damping on the performance of 
miniature mechanical oscillators is a subject of utmost 
importance within the MEMS community from a tech- 
nological standpoint. Additionally, due to the high sur- 
face to volume ratio, flow at the micro/nanoscale shows 
new phenomena not observed at the macroscopic scale 
[21H23| . We carried out a rigorous study of the influence 
of air damping on our devices through a wide range of 
pressures, from 10 mTorr to 760 Torr. 

Two different devices, HI and H2, were studied indi- 
vidually while enclosed in a vacuum chamber, in which 
the pressure was controlled by venting and pumping the 
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FIG. 9: 

Top: Damping coefficient as a function of pressure for device 
HI. Bottom: Damping coefficient as a function of pressure 
for device HI. Calculated curves are shown for the damping 
due to the film (bottom), the top, and the fluid between the 
fingers. 



chamber through a needle valve from 10 mTorr to 1 atm. 
The frequency was swept through the shear mode reso- 
nance after the pressure was stabilized. The pivot mode 
was found to be significantly damped above 40 Torr due 
to the so-called squeeze film damping [ll|. Therefore, 
only the shear mode was studied. The frequency (u)q), 
width (Aw), and amplitude (qo) of the peak were then 
obtained directly from the fit to the curves in Eq. [5] The 
quality factor, Q, was obtained by Q = loq/Auj. The 
frequency and quality factor are plotted in Fig. [8] for de- 
vices HI (top) and H2 (bottom) as a function of pressure. 
The curve for Q is essentially fiat from 10 to 100 mTorr 
since the damping is mainly dominated by intrinsic dis- 
sipation. The quality factor then starts to decrease in 
an almost linear manner until P ^ 100 Torr. At this 
point the quality factor begins to saturate again since in 
the hydrodynamic regime the damping becomes indepen- 
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FIG. 10: 

Damping coefficient as a function of pressure for both 
devices, HI and H2, as a function of K n . The vertical 
dashed lines separate the hydrodynamic, crossover and 
ballistic regimes. 



dent of pressure. However, the resonance frequency stays 
almost constant at pressures below ~ 4 Torr, and then it 
starts to decrease continuously as the pressure increases. 
This behaviour suggests a transition from a highly rar- 
efied (ballistic) regime into a hydrodynamic regime. 

The coefficient of damping, 7, can be extracted ex- 
perimentally from the width of the peak 7 = mAw. The 
intrinsic damping was obtained for both devices from the 
constant value of the width at low pressures. The experi- 
mentally determined damping coefficient as a function of 
pressure is shown in Fig. |5] for HI (top) and H2 (bottom). 
In the two sets of data shown in Fig. [5] the intrinsic damp- 
ing coefficient has been subtracted. The data therefore 
represent exclusively the damping coming from the fluid. 

We followed the recipe given by Bruschi et al. to es- 
timate the theoretical value of the damping 24]. The 
values of the damping coming from different parts, such 
as the center plate and the comb electrodes, can be cal- 
culated using their phenomenological expression for the 
damping of an oscillating plate in close proximity to a 
substrate, the so-called slide film damping: 



l + 2K n 



(8) 



where A p is the area of the movable plate, d is the gap 
between the plate and the substrate, 770 is the viscosity 
at 1 atm, and K n is the Knudsen number, which is the 
ratio of the mean free path to the size of the gap. K n 
quantifies the rarefaction of the fluid, thus the term in the 
parentheses can be understood as an effective viscosity. 
To account for the finite size of the different structures of 
the device, an effective gap size was used. The effective 
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FIG. 11: 

Top: Absorption and dispersion curves for an HI device immersed in liquid 4 He at different temperatures and at 2 bar 
pressure. Bottom: Curves for an H2 device immersed in liquid 3 He at different temperatures and at 3 bar pressure. 



gap is given by \2l 



Resonance as a Function of Temperature in 
Liquid Helium 



l + 8.5d/V 



(9) 



where I is the dimension of the oscillating structure along 
the motion of the plate. The slide him damping coming 
from the hngers (dotted line) and the fluid layer on top 
of the device (dash-dotted line) are almost two orders of 
magnitude smaller than the contribution from the fluid 
trapped between the gap (dashed line); see Fig. HI Thus, 
the total damping coefficient (solid line) is dominated 
almost exclusively by the damping coming from the fluid 
between the plates. These results are consistent with 
those obtained by Bruschi et al. [24j . Furthermore, when 
the damping coefficient is normalized by the geometrical 
constant A p /d, both devices with different gaps exhibit 
a universal behaviour as a function of K n , accross the 
various flow regimes, from hydrodynamic to ballistic (see 

Fig. m 



Low temperature studies in liquid 3 He and in liquid 
4 He were conducted using a dilution refrigerator with 
a base temperature of 5 mK. The chip containing the 
MEMS device was housed in a low temperature exper- 
imental cell made of copper. The cell was loaded with 
the corresponding gas hypercritically to avoid damages 
from capillary condensation. The temperature was var- 
ied by applying heat to the mixing chamber of the di- 
lution refrigerator. At each new heat setting the system 
was allowed to reach thermal equilibrium for a few hours. 
The temperature was measured using a calibrated ruthe- 
nium oxide thermometer located on the mixing chamber 
of the dilution unit. Resonance properties were studied 
in both superfluid 4 He and normal liquid 3 He below 1 K. 
Resonances were obtained for the shear mode of an HI 
device in 4 He while an H2 device was used in 3 He. All 
the measurements at low temperatures were conducted 
in "1/" mode using appropriate levels of dc bias. This 
scheme is advantageous over the "2/" mode in terms of 
heat dissipation, which is a crucial aspect in low tern- 
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perature experiments. The main source of Joule heating 
is from the electrical connections from the bonding pads 
and through the serpentine springs (see Fig. [1]). Using 
the capacitance of the device and the calculated resis- 
tance of electrical connections from the bonding pad to 
the device, the Joule heating in this device is estimated 
to be w 3 fW for 1 V pp ac excitation. We also have con- 
ducted measurements in superfluid 3 He down to 300 fiK 
and no appreciable heating was observed. 

The absorption and dispersion curves are shown in 
Fig. [11] for superfluid 4 He (top panels) between 250 and 
707 mK and for normal liquid 3 He (bottom panels) be- 
tween 56 and 696 mK. A strong temperature dependence 
on the resonance was observed. In the case of 4 He, the 
width of the resonance decreases as the temperature de- 
creases. Below 1 K, the main mechanism of dissipation 
in superfluid 4 He is through phonon excitations, which 
decrease as T 4 in this temperature range. In contrast, 
the width of the peak in 3 He increases as a function of 
temperature. Below its Fermi temperature (1 K), liquid 
3 He begins to display properties of a Fermi system. In 
the Fermi liquid regime, the viscosity of 3 He increases 
very rapidly as ~ 1/T 2 , which explains the broadening 
of the resonance as the temperature decreases. Further 
details of these experiments will be published elsewhere. 
Some preliminary results in superfluid 4 He were discussed 
previously [26j . 

V. CONCLUSION 

A new experimental tool based on MEMS oscillators 
was developed for use in low temperature physics exper- 



iments. This work provides a detailed roadmap of the 
development process from conception, design, and fabri- 
cation to application. Details on the experimental setup 
and the frequency response are presented. The qual- 
ity of the surface was studied in detail by using AFM. 
Two length scales were identified from the autocorrela- 
tion function at 0.137 and 0.591 /zm. From a height his- 
togram the grain size was estimated to be 10.4 nm. 



The devices were studied as a function of pressure be- 
tween 10 mTorr and 1 atm. This allowed us to study its 
resonance properties from a ballistic to a hydrodynamic 
regime. A model of slide film damping was found to be 
in excellent agreement with our experimental data. We 
also demonstrate the use of these MEMS oscillators in a 
cryogenic environment and obtained resonance peaks of 
two devices immersed in liquid 4 He and 3 He below 1 K. 
The devices show great potential in a wide range of ex- 
periments at low temperatures. In particular, they might 
provide a powerful tool to study properties of liquid he- 
lium at the micro/nanoscale. 
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